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Rapidly Rotating Lenses: Repeating features in the 
lightcurves of short period binary microlenses 
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ABSTRACT 

Microlensing is most sensitive to binary lenses with relatively large orbital separations, 
and as such, typical binary microlensing events show little or no orbital motion during 
the event. However, despite the strength of binary microlensing features falling off 
rapidly as the lens separation decreases, we show that it is possible to detect repeating 
features in the lightcurve of binary microlenses that complete several orbits during the 
microlensing event. We investigate the lightcurve features of such Rapidly Rotating 
Lens (RRL) events. We derive analytical limits on the range of parameters where these 
effects are detectable, and confirm these numerically. Using a population synthesis 
Galactic model we estimate the RRL event rate for a ground-based and space-based 
microlensing survey to be 0.32/b and 7.8/b events per year respectively, assuming 
year-round monitoring and where /b is the binary fraction. We detail how RRL event 
parameters can be quickly estimated from their lightcurves, and suggest a method to 
model RRL events using timing measurements of lightcurve features. Modelling RRL 
lightcurves will yield the lens orbital period and possibly measurements of all orbital 
elements including the inclination and eccentricity. Measurement of the period from 
the lightcurve allows a mass-distance relation to be defined, which when combined 
with a measurement of microlens parallax or finite-source effects, can yield a mass 
measurement to a two-fold degeneracy. With sub-percent accuracy photometry it is 
possible to detect planetary companions, but the likelihood of this is very small. 

Key words: gravitational lensing: micro - celestial mechanics - binaries: general - 
planetary systems - Galaxy: bulge 



1 INTRODUCTION 

By monitoring hundreds of millions of stars towards the 
Galactic bulge and Magellanic cloud s. Gravitationa l mi- 
crolens ing surveys such as OGLE (|Udalskil 120031 ) and 
MOA (|Hearnshaw et aLllJOoil ) detect ~ 1000 microlensing 
events per year. T he lightcurv e s of m ost microlensing events 
follow the typical IPaczvriskil (|l986l ) form of a point-mass 
lens with a point source. However, many event lightcurves 
have a more complex form due to th e effects of a bi- 
nary o r planetary companio n to the lens (jMao fc Paczvriskil 
I1991I : iGould fc Loebl Il992). a binary companion to the 
sourc e (iGriest fc Hul Il9921 ). microlens par allax ( Refsdal 
19661: [Gouldl Il992l). finite-source size llGouldl Il994 



Witt fc Maolll994l : iNemiroff fc Wickramasinghdll994l ). or a 
combination of such. These effects provide additional in- 
formation about the lens, such as the mass ratio and pro- 
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jected separation in binary or planetary lens events, or mass- 
distance relationships with parallax or finite-source effects. 
Measurement of both the finite-source size and microlensing 
paral lax allow the lens mass to be solved for uniquely d Go uldl 
Il992h . and if these occur in a binary or planetary lensing 
event, also allow measurement of the companion mass. 

The complexities of microlensing lightcurves, to a 
greater or lesser extent, can all be considered as devia- 
tions from the single lens Paczyhski form. The deviations 
may be relatively minor and can co ver the entire lightcu rve, 
as in most parallax events (e.g. ISmith et al.l l2002al ). or 
they can be large and cover only a small fr action of the 
lightcurve, as in many bin ary lens events (e.g. iKubas et al.l 
I2OO5I : iBeaulieu et al.l 120061 ). In binary lens events, these de- 
viations from the single lens form are caused by a difference 
in the magnification pattern of the lens. The most promi- 
nent 'features' of the binary lens magnification pattern are 
caustics, where the magnification of a point source diverges 
(see Figure [1]). A source passing over a caustic will show a 
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sharp rise in its magnification as it enters it and a sharp 
fall as it leaves. Other, more smooth magnification pattern 
'features' can be associated with the caustics in some way. 
The caustic features of a binary lens magnification pat- 
tern are a natural feature of the binary lens mapping 
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which maps the angular position of the source to image posi- 
tions under the influence of the lens, a nd where w e have used 
complex coordinates (e.g. z = x -\- it/. IWittlll99d 'l: bars rep- 
resent complex conjugation, Zb is the position of the source, 
z the position of the image, z\ and Z2 the positions of the 
primary and secondary lens components respectively, and 
q = M2/M1 is the mass ratio of the lens components. All 
angles have been normalized to the Einstein ring radius 



TE 1 /4G ,, N r^ ,^ 
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where te is the physical Einstein radius, Di and D^ the 
distance to the lens and source respectively, M — M1 + M2 is 
the total lens mass, x — Di/Ds is the fractional lens distance, 
and G and c are the gravitational constant and speed of light 
respectively. The magnification fi of an image is given by the 
determinant of the Jacobian of the lens mapping 

(3) 
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The magnification diverges when det J = 0, and the so- 
lutions of this equation form smooth, closed curves in the 
image plane called critical curves, which when mapped back 
to the source plane form closed, cuspy curves: the caustics 
(see Figures [T] and (Sjl . 

In a binary lens event, the caustics are largest when the 
projected lens separation s — \z2 — zi\ ~ 1, i.e., the lens com- 
ponents orbit with semimajor axis a ~ te ~ 2-3 an. At these 
separations there is only a single, so-called resonant, caus- 
tic. Lenses with wider separations have two caustics, while 
those with closer separations have three ({Schneider fc Weisa 
119861 ). The caustics of close and wide systems have smaller 
caustics than resonant systems, and so the source is less 
likely to encounter them. Far from the caustics, the binary 
lens features tend to be weak, and therefore, the lightcurves 
of binary lenses with very close or very wide orbits (which 
have correspondingly very small caustics), in most cases, 
will be i ndistinguishable from those of a single lens (e.g. 
iGaudi fc Gould 1997). 

The binary lenses with the strongest lightcurve features, 
thus have orbital periods T ~ 1000 d, much longer than the 
microlensing event timescale 
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for a typical Galactic microlensing event, where vt is the rel- 
ative lens-source transverse velocity (source velocity). The 
lenses therefore complete only a small fraction of their orbit 
during the course of the microlensing event, and only a frac- 
tion of the events are expected to s how detectable signs of 
orbital m otion i n their lightcurve s (IGaudi fc Gouldl 1997 : 



pominil4_ 19981: Konno fc Koiima 19991: loka et al 



199c 



Rattenburv et al.l I2OO2I : IPennv et all bOllT ). Those events 



where it has been detected typically involve caustic cross- 
ings, because the sharp caustic crossing features on the 



lightcurve allow precise timin g, which can be used to con- 
strain even small lens motions ( Albrow et al. 2000l: An et al.l 
'2002'; 'Gaudi ct al."2008'; 'Rvu ct al."2010 HHwang et al.ll2010r 
Skowron et al. 2011; Batista ct al. 201l|). In these events, 
the orbital motion has constrained the properties of the lens 
orbit, and in one c ase allowed the mea surement of several 
orbital parameters (jBennett et al.ll201(]r ). 

In a previous paper, we have shown that with accu- 
rate photometry and dense enough lightcurve coverage, it 
is possible to detect orbital motion in lenses w ith closer or- 
bits, and weaker, smooth lightcurve features (jPennv et al.l 
l2011f ). Typically occurring in events with longer microlens- 
ing timescales, t-^ ~ 100 d, the shorter orbital periods means 
that the lens completes a much larger fraction of the orbit 
during the event. In this paper we investigate the extreme 
case, where the lens completes at least one orbit during the 
microlensing event. In events involving such rapidly rotat- 
ing lenses, magnification pattern features can sweep over the 
source more than once, and if detected in the lightcurve, 
these repeating features allow an accurate measurement of 
the lens period. Knowledge of the period places constraints 
on the lens mass, and if combined with a measurement of 
the finite-source effect or parallax, can be used to mea sure 
the lens mass to a two- fold degeneracy (|Dominikll 19981 ) . 

We begin the paper by defining when a binary lens is 
a rapidly rotating lens in Section (2] In Section [3] we look 
at whether rapidly rotating lenses are detectable, and the 
rate at which they are expected to occur. In Section UJ 
we then look at what information can be gained from ob- 
serving a rapidly rotating lens, and suggest a method for 
modelling such events. In two appendices we examine how 
rapidly rotating lenses affect the microlensed images, and 
introduce additional effects that can affect rapidly rotating 
lens lightcurves. 



2 WHAT IS A RAPIDLY ROTATING LENS? 

We define a rapidly rotating lens (RRL) to be a binary mi- 
crolens, which if monitored continuously with suitable pho- 
tometric accuracy, would guarantee that at least one feature 
of its magnification pattern would be seen to repeat at least 
once in its lightcurve, due to the lens orbital motion. This 
implies that the lens completes at least two orbits during 
the time in which its binary lensing features are detectable. 
We choose this definition over the more simple comparison 
of microlensing and orbital timescales (e.g. T < Ie) because 
without detecting binary features it is impossible to measure 
the binary's rotation. As mentioned in the previous section, 
the strength of binary features declines as the orbital sepa- 
ration and period decrease, so simply decreasing the period 
does not necessarily increase the prospects of detecting a re- 
peated feature. A rapidly rotating lens must therefore com- 
promise between a fast rotation rate and detectable binary 
lensing features. 

Throughout the paper we shall focu s on close topology 
lense s , which h ave separations s < 0.7 ({Schneider fc WeissI 
119861 : lErdl fc S chneider 1993), a choice we shall justify in 
Section [3] Figure [1] shows the magnification pattern of a 
close topology lens, and labels a number of features. The 
structure and features of the magnification pattern depend 
only on the projected separation of the lens components s. 
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Figure 1. The magnification pattern of a close topology mi- 
crolens. The dots denote the lens positions, with the primary lens 
at negative x. The lens has a mass ratio q = 0.3 and projected 
separation s = 0.6. Notable features of the magnification pattern 
are labelled. 
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Figure 2. The lightcurve of a rapidly rotating lens. The up- 
per panel shows the RRL lightcurve in black, and the Paczynski 
lightcurve of a single lens with the same total mass in grey. The 
lower panel shows the residual with respect to the Paczynski 
lens lightcurve. Features due to the magnification arms appear 
as peaks in the residual, while between them there are relative 
demagnifications. Large, short duration spikes occur when the 
secondary caustic passes close to or over the source. The system 
has parameters is = 61 d, T = 92 d, s = 0.23, q = 0.8, mo = 0.3, 
<pO = 1.75 (see Section |4] for definitions of «o and (po). 



and the mass ratio q (JErdl fc Schneider! [l993J) . The most 
important features of the close magnification pattern are 
a central caustic, located at the lens centre of mass, and 
two secondary caustics which lie away from the lens centre. 
Stretched between the central and secondary caustics are 
two 'arms' of excess magnification (relative to the magnifi- 
cation that would be caused by a single lens of mass equal 
to the total binary mass). During a microlensing event, a 
source will travel across the magnification pattern, and we 
will observe the source to change in brightness, the form of 
this lightcurve being determined by the trajectory that the 
source takes. As the source moves, the magnification pattern 
will not stay fixed, as the binary will also move in its orbit. 
Should the lens orbit lie face-on to the line of sight, then the 
magnification pattern will rotate as the source moves across 
it. Should the orbit be inclined or eccentric, the structure of 
the magnification pattern will a lso change, as it depends on 
the projected lens separation s ([Schneider fc WeisJI 198(3 ). 

Figure[2]shows the lightcurve (thick black line) of a RRL 
with a similar magnification pattern to that shown in Fig- 
ure [T] It closely resembles the lightcurve of a single lens, the 
Paczyriski lightcurve (grey line), but with a quasi-periodic 
variation over the entire lightcurve that only becomes ob- 
vious in the residual that is left once the Paczynski curve 
is subtracted from the lightcurve. These periodic features 
correspond to the magnification arms that extend between 
the secondary and central caustics, which sweep over the 
source as the lens rotates. The microlensing timescale of the 
lightcurve shown is Ie = 60 d, but it is clear that repeating 
binary features remain in the lightcurve at time from peak 
magnification much greater than this, which corresponds to 
a source position far outside the Einstein ring. This is be- 
cause the secondary caustics can lie far outside the Einstein 



ring, their distance from the lens centre increasing as the 
binary separation decreases. However, their size, and the 
strength of the magnification arm connecting them with the 
central caustic, decrease with decreasing binary separation. 
We note at this point, that despite the large separation of 
the secondary caustics, we nee d not consider relativ istic ef- 
fects of superluminal caustics (jZheng fc Gouldll200(]l l as the 
ratio of the caustic rotational speed to the speed of light in 
all the cases we will consider is ~ 10~^. 

An RRL can clearly exhibit interesting, repeating 
lightcurve features if the binary period and separation con- 
spire, but until now we have only considered RRLs in re- 
lation to parameters normalized to the typical microlensing 
length- and time-scales. In order to see if RRL events will 
be detectable in real microlensing surveys we must consider 
how their properties relate to the physical parameters of the 
lensing system. 



3 ARE RRLS DETECTABLE? 

In the previous section we defined a criterion for a lens to 
be a rapidly rotating lens, and described the features of an 
RRL event. In this section we put the definition on a more 
quantitative basis, and investigate whether RRLs will occur 
amongst the microlenses that are probed by microlensing 
surveys. To determine if detection is plausible, we investi- 
gate the range of physical parameters required to produce 
a microlensing event with repeating features, first analyti- 
cally and then numerically. Finally we apply our numerical 
method to simulated microlensing surveys to estimate the 
expected rate of RRL detections. 
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3.1 An analytical approach 

To see repeating features in a microlensing event, tiie most 
fundamental requirement of the system is that the lens com- 
pletes more than one orbit during the event. The magnifica- 
tion pattern of a binary lens is complicated, but the essential 
features of a close binary lens can be captured by assuming 
it to be composed of two straight, radial arms that extend 
from the centre of mass to the position of the secondary 
caustics. Under this assumption, and assuming a random 
initial phase angle, repeating features are guaranteed to be 
observed if the lens completes two orbits in the time that 
the source spends within the radius swept out by the arms. 
We can write this as an inequality 



2T < -u±tE, 



(5) 



where u± is the radial position of the secondary caustics in 
units of Einstein radii (see Figure [8]), and the factor of j is 
the mean chord length across a unit circle, and accounts for 
the random impact parameter of source trajectories relative 
to the lens centre of mass. It should be noted that it is pos- 
sible for a feature to repeat if the binary completes between 
one and two orbits, but this requires a coincidence in the 
timing of the first feature. 

Both the orbital period and the Einstein timescale de- 
pend on the lens mass, and the period also depends on the 
lens semimajor axis, so it is possible to write this constraint 
in terms of M and a. For projected lens separations s <^ 1, 
iBozzal (|2000|) has derived an analytical approximation for 
the secondary caustic positions (see Equation I18|l. which if 
we keep only the first order terms is 



u±{s,q) 



(6) 



Using equations [2] and |4] and Kepler's third law, with a little 
algebra we can then write Equation [5] as a constraint on the 
semimajor axis of the binary 



a < 4.51 AU [x{l - x) 
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where we have assumed a face-on orbit so that s = a/r-E, 
and where M is the total lens mass in solar masses, D^ the 
source distance in kpc, x = Di/Ds is the ratio of lens and 
source distances and vt the relative lens-source velocity in 
km s~^. 

While we have an upper limit on the lens semimajor 
axis, in order to detect the RRL lightcurve features, they 
must be strong enough to be detectable in the photometry 
of the microlensing event. This requirement is somewhat am- 
biguous, but as the magnification pattern depends only on s 
and q, and the strength of features decreases with decreasing 
s, we can assume that for a given photometric precision and 
mass ratio, magnification pattern features will be detectable 
only when the separation is larger than a certain value, i.e. 



S > Sdctj 



(8) 



where Sdet depends on q and the photometric accuracy. For 
stellar binary mass ratios, there will only be a small depen- 
dence on q, but there will be a strong dependence on the 
photometric accuracy; a value of Sdct = 0.3 is reasonable, 
however (see Section[T2J ■ We can again write this constraint 
as a limit on the semimajor axis 




a > 2.85 AU Sdct[2;(l - x)]^^'^dI^^ M^^^ . 



(9) 
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Figure 3. Plot showing the region of the total mass-semimajor 
axis plane where repeating features are observable. The red line 
shows the upper limit on a, provided by the constraint in Equa- 
tion [T] while blue lines show the lower limit on a, provided by the 
constraint in Equation |9] with values of Sdet = 0.3,0.2 and 0.1 
from top to bottom. The other parameters are set at x = 0.75, 
Db = 8 kpc and vt = 50 km s~^. The region where repeating 
features are detectable for Sdet = 0.3 is shaded grey. 



We now have two constraints on a, an upper and a lower 
limit, which are dependent on other parameters of the lens- 
ing system, the most interesting being the total lens mass. 
Figure |3] shows the two constraints on semimajor axis as a 
function of mass, for a lens system with x = 0.75, D^ = 8 kpc 
and Vt — 50 km s^^, with values of Sdet = 0.3,0.2 and 0.1. 
Other than the slow lens-source velocity {{vt) ~ 200 km s~^ 
for a Bulge microlensing event), these values are typical of 
a microlensing event towards the Galactic Bulge. The plot 
shows that most of the parameter space is excluded, but 
thanks to the differing power-law indices on the mass de- 
pendence, there is a small range of parameters over which 
repeating features should be detectable. For the parameters 
shown, the 'detectable region' opens up at M ~ ^Mq and 
a ~ 1 au, and widens to 3.3 < a/au < 4.4 by M = IOMq. 
The dependence of the limits on other parameters means 
that the region of detectability will get smaller, and move 
to larger a as the source distance grows, will get larger and 
move to smaller a and M as the lens moves closer to the 
source or the observer, and will get smaller as the relative 
lens-source velocity increases. A small but significant frac- 
tion of binary stars will have total masses and semimajor 
axes in the range of detectability (e.g. iDuguennov fc Mavorl 
Il991f ). especially if improved photometric accuracy can re- 
duce Sdct- 



3.2 A numerical approach 

In deriving analytical limits on the range of lens parameters 
we have had to make assumptions about the magnification 
pattern and strength of features. If we instead proceed nu- 
merically, we need not make these assumptions as we can 
determine precisely the regions of the magnification pattern 
where features are detectable for any given photometric ac- 
curacy. We define a detectability e that is the probability 
that, for a given lens system and photometric precision, an 



Rapidly Rotating Lenses 5 




M/Mq 

Figure 4. Plot of the average detect ability e against total lens 
mass M and seniimajor axis a, for a lens with mass ratio q = 
0.3. The lens and source distances and relative velocity are the 
same as in Figure \3\ The red line is the analytical upper limit 
of Equation [7] while the blue line is the analytical lower limit of 
Equation [9] with a value of Sjet = 0.28 for a photometric precision 
^"^min = 0.01. The green line at the top of the figure marks the 
boundary between regions of close and resonant topology lenses 
— we only calculate e for close topology lenses. The green points 
lower in the figure mark points where our calculation of e failed 
(see text for details). 



RRL with a face-on orbit will exhibit at least one detectable 
repeating feature in its lightcurve. A feature is said to be de- 
tectable at a radial position u, if the range of magnifications 
fj. over a circle of radius u satisfies 



Am = 2.5 log 



^J■n 



Mmin(w) 



> Arun 



(10) 



where we have expressed the range of magnifications 
(Mmin — >■ Mmax) as a magnitude difference Am, and Arrimin is 
the photometric detection threshold, which can be taken to 
mean the typical uncertainty in magnitude on a data point 
in the baseline of the lightcurve. In calculating e we average 
over the random parameters of the source trajectory and 
phase angle. 

We can now test the predictions we made in Section FS.ll 
by mapping the detectability e against total mass M and 
semimajor axis a, for the set of parameters we used for Fig- 
ureO Figure |4] shows regions of finite detectability in shades 
of gray up to black when e = 1. Green points in the plot show 
where the calculation of the detectability, which requires sev- 
eral numerical minimization and root finding steps, failed. 
The green line at the top shows the boundary between the 
close and resonant caustic regimes, where the three caustics 
of the close topology merge into a single resonant caustic. 
We do not calculate the detectability in the resonant regime. 
The red and blue lines show the analytical upper and lower 
limits of Equations [7] and [H however with Sdot = 0.28 as 
opposed to 0.3. It can be seen in the figure that the ana- 
lytical upper limit of Equation [7] agrees very well with the 
numerical region of detectability, coinciding with the bound- 
ary where e begins to fall from unity with increasing a over 
the entire range of M shown. Equation[7l without the factor 
of 2 that was introduced on the left-hand side of Equation [5] 
to guarantee a repeated feature, also describes well the re- 



gion where detection becomes possible but is not guaranteed 
(i.e. <£< 1). 

The analytical lower limit, once the parameter Sdct has 
been adjusted to 0.28 for a guaranteed repeating feature, 
also agrees well with the numerical region of detectability. 
It should be noted however, that the slope of the lower edge 
of the numerical region is slightly shallower than the analyt- 
ical lower limit. This becomes more pronounced when the 
lens gets closer to the source, the total mass increases or 
the source velocity decreases. This is because the assump- 
tion, that there are detectable features over the entire mag- 
nification pattern within u < w±, breaks down, and the 
detectable features lie in two disjoint regions, a disc sur- 
rounding the central caustic and an annulus containing the 
secondary caustics. The size of these regions depends on s, 
and the lower limit on a (Equation |9]) becomes a shallower 
function of A/. This effect is more important in determining 
the slope of the lower limit on a where e = 0. 

Having looked at the role of mass and orbital separa- 
tion, it is important to investigate how the detectability of 
repeating features depends on other factors. Figure [5] shows 
detectability maps similar to that in Figure lU but for a mass 
ratio q — 0.99 and various values of the lens distance, source 
velocity and photometric detection threshold. The maps are 
arranged into four panels with different photometric thresh- 
old values of Ammin = 0.005, 0.01, 0.02 and 0.04 from left 
to right, top to bottom. Each panel is made up of a grid of 
a-M maps showing detectability for different fractional lens 
distances x in each column, and different source velocities 
in each row. Looking first at the grid with Ammin ~ 0.01 
(top right), it is clear that the source velocity has a large 
effect on the detectability, with large regions of detectability 
for vt = 50 km s~^ at all lens positions, which are reduced 
drastically for Vt = 100 km s~^. Once vt ~ 150 km s~^ 
there is only a tiny region of low detectability for lenses 
very close to the source (or to the observer, as x{l — x) is 
symmetric about x = 0.5). For vt — 200 km s~^ there is only 
detectability in the most favourable cases of very high pho- 
tometric accuracy and fractional lens distance. This strong 
dependence on Vt occurs because the number of orbits com- 
pleted by the lens decreases as Vt increases (the v^ term 
in Equation [7| but does not affect the strength of binary 
features (Equation |9] is independent oi Vt) - there is no can- 
cellation of the Vt term in the ratio of upper to lower limits, 
while there is a partial cancellation for all the other phys- 
ical parameters. Unfortunately, the microlensing event rate 
peaks at vt ~ 200 km s~^, but there is a significant fraction 
of events with vt < 100 km s"^ (e.g. iDominikI [20061 '). 

The lens distance does not affect the size of the de- 
tectable region as strongly as the source velocity does, as 
the upper and lower limits of the detectable region scale with 
x{l — x) as similar power laws (—0.4 and —0.5 respectively). 
However, this similar scaling means that the detectable re- 
gions move as x changes, occurring at lower a, and increasing 
in size slightly, as the product x{l — x) decreases. For mi- 
crolensing events to wards the Galac tic Bulge, the event rate 
peaks at x ~ 0.8 (e.g |Dominikll2006l l . whereas for self lensing 
in the Magellanic clouds x will be close to unity, x ~ 0.98. 

The photometric precision of the observations strongly 
affects the detectability of repeating features. For Ammin = 
0.005 and 0.01 we see large regions of detectability for small 
source velocities, and for Ammin ~ 0.005 even some de- 
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Figure 5. Maps of detectability £ plotted against a and M for a binary of mass ratio q = 0.99 and various values of the fractional lens 
distance x = Di/Ds, source velocity vt and photometric detection threshold Ammin. Each small sub-panel is essentially the same as the 
plot in Figure[4l but with different parameter values and a slightly restricted range 0.1 ^ M/Mq ^ 3, 0.1 ^ a/au ^ 3. From left to right, 
top to bottom panel have increasing values of the photometric threshold Ammin = 0.005, 0.01, 0.02 and 0.04. Moving from left to right, 
sub-panels have different fractional lens distances x = 0.5, 0.75, 0.9, 0.95 and 0.98; the results remain the same under the transformation 
X —^ (1 — x), i.e. there is reflectional symmetry about x = 0.5. Moving from bottom to top, sub-panels have different source velocity 
Vt = 50, 100, 150 and 200 km s~^. The source distance is fixed at Db = 8 kpc. The black line shows the boundary between close and 
resonant caustic structures, above which we do not plot e. As in Figure |4] there are points where the calculation of e fails, but these are 
not shown for clarity, as they do not impinge on the regions of detectability. 



tectability when Vt = 200 km s^^. As Arnniin increases to 
0.02, the detectable regions shrink significantly and all but 
disappeared for Vt ^ 100 km s~^. For Ammin = 0.04 there 
is virtually no detectability, with only a small chance of de- 
tection for the smallest velocities and largest lens distances. 
Increasing the threshold effectively increases the lower limit 
of a at which binary features are detectable, while leav- 
ing the upper limit unchanged, and like the source veloc- 
ity, the photometric threshold has a large effect on the size 
of the detectability region. It should be noted that the de- 
tection threshold Ammin is in fact a combination of the ef- 
fects of photometric precision and the blending by unrelated 
starlight, which acts to add a noise component to the mea- 
surement of the magnification caused by the lens. The effect 
of blending is discussed further in Appendix [Bl] 

Even in the most favourable case of low photometric 
threshold and source velocity, and high fractional lens dis- 
tance, the region of detectability does not reach the bound- 
ary between close and resonant caustic topologies. This 
is because as the projected separation increases, near the 
topology boundary, the secondary caustics move inwards 



rapidly before merging with the central caustic, decreasing 
the range over which binary features are detectable. At the 
same time the orbital period will increase rapidly as the 
semimajor axis increases. These combined effects mean that 
in order to see repeating features from a lens with resonant 
topology, an extremely low source velocity is necessary to 
allow the lens to orbit in the time the source spends near 
the resonant caustic. 

Figure|6]shows the same maps as Figure[5]but for differ- 
ing q, and the threshold fixed at Ammin ~ 0.01. The maps 
for q = 0.3 are similar to those for q = 0.99, and there 
is little difference in the size of the region of detectability. 
However, once q has fallen to 0.1, the size of the detectable 
region has begun to shrink, such that for higher values of 
Ammin (not shown) there is only a very small chance of de- 
tection with small source velocities. For lower mass ratios 
still, there are only very small regions of detectability for 
g = 0.03 and effectively zero detectability for g = 0.01. If we 
take the boundary between brown dwarfs and planets to be 
at 13Mjupiter, there is a very small region of detectability 
where the secondary lens is a planet, but the point of the 
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q = 0.3, A'm„i„ = 0.01 

X 0.5 0.7.5 




q = 0.03, Am„,i„ = 0.01 

X 0.5 0.75 



Figure 6. As Figure[5l but for differing mass ratios. Moving from top left to bottom right the detectabihty is plotted for q = 0.3, 0.1, 0.03 
and 0.01. The total mass corresponding to a secondary below the deuterium burning limit, Myi ~ 13A^jupiter, is M < 0.054Mq, 
M < 0.14AfQ, M < 0.43Mq and M < 1.25M0 respectively for each value of q. The photometric detection threshold in each case is 
Am^in = 0.01. 



detectable region where the upper and lower limits meet oc- 
curs close to this point, regardless of the mass ratio. So, there 
is little chance of detecting repeating features from a plane- 
tary system, unless the photometry is very accurate, the lens 
very close to the source, or the source velocity is significantly 
smaller than 50 km s~^. Such low velocity events are rare, 
but are known to occur, e.g. the event OGLE-1 999-BLG-19 
had a source velocity vt = 12.5 ± 1.1 km s~^ (jSmith et al.l 
l2002bl ). 



3.3 How many RRL events will we detect? 

To estimate the RRL event rate we conducted a simula- 
tion of a space-based //-band m icrolensing survey, such 
as E uchd ^caulicu et al.1 |2010| ) or WFIRST (Bennett 
[20l2), and a ground-bas e d J-b a nd surv e y, ba sed on 
OGLE-III l|Udalski et al.l 1 19971 : lUdalskil |2003| ). Us- 
ing the Besangon popul ation synthesis model of the 
Galaxy (| Robin et al. 2003) , including a three dimensional 
extinction model (Marsh all et al.l l2006r ). we produced a 
catalogu e of possible micr olensing events following the 
recipe of lKerins et al.l (|2009l ). Source stars with magnitudes 
Hs < 25 and 7s < 21 are drawn from the Besangon model 
and lensed by stars of any magnitude in the space- and 



ground-based simulations respectively. The lens mass is 
split up into two components with a mass ratio q distributed 
logarithmically in the range 0.1 ^ q < 1, and orbit with 
a semimajor axis a distributed logarithmically in the 
range 0.1 ^ a/au < 4. Each event is assigned a weighting 
w = 2rBi'tWo,max proportioual to its event rate, where 
uo,max is the maximum impact parameter that the event 
could have and its peak single lens magnification remain 
detected at 5a" above baseline, taking into account blending. 
Each event was assigned a blending f raction fg ^ 1 draw n 
from the blending distributions of ISmith et al.l ((20071), 
with source density 131 stars per square arcmin, and 
seeing 0.7 arcsec and 1.05 arcsec for the space-based and 
ground-based simulations respectively. This will overesti- 
mate the blending effect for the space-based simulation, as 
the diffraction limited PSF for aim telesco pe will have 
a full width half maximum ~ 0.4 arcsec, but ISmith et al.l 
([2007J) do not simulate seeing better than 0.7 arcsec. The 
final blending suffered by the source /a also includes fiux 
from the lens, which is obtained from the Besangon model 
assuming it is a single star. The severity of blending is 
thus overestimated, especially for the space-based survey, 
and as blending has a large effect on the detectabihty (see 
Appendix IB1[) . the event rates we estimate will be con- 
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servative. However, we do not include the efTect of orbital 
inclination, which can decrease the amplitude of lightcurve 
features slightly (see Appendix IB3[l . so this optimistic 
assumption will likely balance the pessimistic blending we 
apply. The photometric detection threshold was calculated, 
based a pproximately on the proposed desi gn of the Euclid 
mission (jEuclid Pavload Definition Docume nt 2009 .) for the 
space -based survey, and the OGLE-HI setup (|Udalski et alj 
[l997|) for the ground-based setup. Total event rates are 
normalized to rates Ful = 7000 yr~^ for the space-based 
survey fe.g. ISennett fc RhidlJOO^ '). and T^^l = 600 yr"^ for 
the ground based survey, corresponding roughly with the 
OGLE-HI survey. The rate of RRL events Frrl is taken to 
be 



Ff 



FpL 






(11) 



the normalized sum of the product of Zi and Wi, the de- 
tectability and weight of event i respectively, over all mi- 
crolensing events, where W = ^}^Wi again summed over all 
events. 

The simulations do not account for the observing strat- 
egy, and assume that frequent monitoring (a few data points 
per night or greater) is conducted for a significant fraction 
of the year (6 months or greater) . It is difficult to assess the 
impact of seasonal observability on the probability of de- 
tecting repeating features, without performing detailed de- 
tection efficiency simulations. To account for this we intro- 
duce a factor /soas , the fraction of a year spent continuously 
observing, which is approximately the probability that an 
individual feature is 'caught'. We m ust also account f or the 
fact that not every lens is binary. iRaghavan et al.l (|20ld ) 
find that 44 percent of stellar systems are multiple, with 
mass ratios q > 0.1, and of these about 20 percent lie in the 
appropriate semi-major axis range, so we adopt a binary 
fraction /b ~ 0.1. 

For our entire sample of space-based survey events we 
find that RRL events make up a fraction (1.1 ± 0.2) x 10~^ 
of the total microlensing event rate, which corresponds 
to an event rate Frrl = (7.8 ± 1.5)/soas/b yr"^. Simi- 
larly for the ground-based survey we find that a fraction 
(0.5 ±0.1) X 10~^ of the total microlensing event rate is 
made up of RRLs, which corresponds to an event rate 
Frrl = (0.32 ± 0.06)/scas/b yr~^. In all cases the errors 
are statistical. 

Figure [7] shows the distribution of microlensing 
timescales for the detectable RRL events and all microlens- 
ing events in the space-based simulation. The results are 
very similar for the ground-based survey, other than the 
overall normalization. The distributions do not take into 
account any timescale dependence on detection efficiency, 
or the binary fraction. The timescale distribution for RRLs 
shows a peak at fE ~ 200 d that is at timescales a factor 
of ten longer than the overall microlensing timescale dis- 
tribution. Even at this timescale however, detectable RRL 
events make up less than one percent of the whole. As the 
timescale increases, the fraction of RRL events increases. 
Long timescale events are intrinsically rare, but RRL events 
make up a significant fraction of all events with these 
timescales, and so such events are good targets to search 
for RRL signals. Additionally, their long timescales mean 
that each event is observable for many years, and it is pos- 
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Figure 7. Microlensing timescale distributions for detectable 
RRL events (solid line) and all microlensing events (dashed line) 
for the space-based survey. Other than the overall normalization, 
the distribution for a ground based survey is very similar. The dot- 
dashed lines show the expected asympto tic slope of the timescale 
distribution, with power law indices ±3 IIMao fc Pac2viiskilll996n 



sible to obtain dense coverage of the lightcurve with stan- 
dard survey-mode observations. The timescale distribution 
for all ev ents agrees well with th e expected asymptotic be- 
haviour IjMao fc Paczviiskil 119961 ). except for the points at 
very small and large ts, where small number statistics are 
in effect. 

Various microlensing surveys have targeted the Galactic 
bulge, pretty much continuously for roughly twenty years. 
These survey-mode observations take place over much of 
the year, so the seasonal observability factor /seas will be 
close to unity. There is therefore a good chance that there 
is of the order of one RRL event in current microlensing 
data sets. New ground based microlensing surveys, already 
in operation and due to start in the near future will increase 
the overall microlensing event rate significantly, and so there 
is also a reasonable chance of detecting of the order of one 
RRL over a timescale ~ 5 years. 

A space-based microlensing survey is proposed for two 
space missions which would launch at th e end of the 
decade: E SA's Eudid llBeauheu et al.l [2OI0I ) and NASA's 
WFIRST (lBennettll201ll ). Such a mission may only spend 
2-3 months per year performing a microlensing survey, as 
the majority of observing time would be spent on dark en- 
ergy surveys. As such the seasonal observability /soas ~ 0.2 
factor would be low, as a high degree of coincidence would 
be necessary for multiple RRL features to fall within the 
observing windows. The number of space based RRL detec- 
tions is therefore likely to be low in reality. However, a dedi- 
cated space-based microlensing survey, possibly as a mission 
extension to Euclid or WFIRST, observing the bulge con- 
tinuously for most of the year would be very likely to detect 
RRL events. 



4 PHYSICAL PARAMETERS FROM RRLS 

The lightcurve of a static binary microlensing event contains 
information on the lens, which can be found by fitting the 
lightcurve with a static binary microlensing model. Simi- 
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Figure 8. Parametrization of the rapidly rotating lenses. Caus- 
tics are shown as solid red lines, the lens positions as red circles, 
with the primary lens in the positive quadrant, and the source 
trajectory as a solid green line at positive x. The green line at neg- 
ative X shows the trajectory of a source with negative mq (see text 
for more details). The binary axis (BA), which subtends an angle 
(/>(t) relative to the fixed x-axis, rotates at a frequency u) = 2n/T. 
(u-l-,i/'+) is the position of the secondary caustic in polar coor- 
dinates fixed to the binary axis. The blue dotted line shows the 
iBozzal pOOOh approximation to the position of the centre of the 
secondary caustic (Equation II8I I for this lens. The lens has the 
parameters s = 0.65 and q = 0.1, and lengths are in units of 
Einstein radii. 



larly, the lightcurve of a rapidly rotating lens contains infor- 
mation about the lens and its orbit. In this section we in- 
vestigate the information it is possible to extract from RRL 
lightcurves, and how this can be done. 

The static binary lens lightcurve for a point source, can 
be described with a minimum of seven parameters: three to 
describe the source trajectory, usually an impact parameter 
Wo and angle a, and the time of closest approach to the origin 
to; one for the lightcurve baseline mb; two to describe the 
lens, the mass ratio q and projected separation s in units of 
Einstein radii; and finally the Einstein radius crossing time 
tE. The coordinate system is usually chosen so that both 
lenses lie on the 3;-axis, and the origin is the centre of mass; 
we shall refer to it as the static centre of mass system. 

The simplest RRL, with a face-on, circular orbit re- 
quires just one additional parameter, the orbital period T, 
for a total of eight parameters. In contrast, a full Keple- 
rian orbit requires five additional parameters (including the 
period), bringing the total to thirteen parameters, many of 
which will be hard to constrain. We demonstrate below that 
the eight parameters of the face-on, circular model can be 
well constrained by the lightcurve, and parameters can effec- 
tively be 'read off' the lightcurve with only a small amount 
of algebraic manipulation. It should be possible to use these 
parameter estimates in a more detailed modelling analysis, 
either using the face-on, circular model (which will be well 
constrained, should the face-on, circular orbit approxima- 
tion apply), or as partial constraints for a full Keplerian 
model. This analysis, which we describe briefly later, can 



significantly reduce the range of parameters it is necessary 
to search in order to find a suitable event model. In Ap- 
pendix IB3I we briefly discuss the effects of orbital inclina- 
tion and eccentricity on the lightcurves and detectability of 
RRLs, and in Appendix IB4I we discuss the effect of parallax 
on an RRL lightcurve. 

We choose a coordinate system flxed with respect to 
the sky, with its origin the lens centre of mass. As such, the 
lens components are not flxed. For convenience, we recast 
the angle a — >■ (/f>o, where cj>o is the angle subtended by the 
primary mass relative to the a;-axis at time to, and we flx the 
angle of the source trajectory, such that the source travels 
parallel to the y-axis. At time t the source is at the (complex) 
position 



Zsit) = I UQ 



t-to 

tE 



and subtends the angle 



9{t) — arctan 



to 



UotE 



(12) 



(13) 



with respect to the a;-axis. Similarly, the binary axis, which 
we define as the line extending from the centre of mass 
through the primary mass, subtends an angle 



0W = f(t- 



to) + 00 



(14) 



with respect to the a;-axis. This parametrization is shown 
in Figure m The parametrization differs from that recently 
proposed bv lSkowxon et al.l (|201l|) for orbiting binary lenses, 
which is best suited for binaries wit h orbits much longer 
than the microlensing timescale. The lSkowron et al.l (120111 ) 
parametrization is expressed in terms of the 3-dimensional 
position and velocity of one lens component, as the "on- 
sky" position components will be well constrained, the "on- 
sky" velocity components may be well constrained and the 
radial position and velocity are likely to be poorly or not 
constrained. However, as we will show, for an RRL it is the 
orbital period and phase angles that will be well constrained, 
so it is better to couch the problem in terms of these quan- 
tities. 

Many of the features in a close binary lens magnifica- 
tion map are radial, or approximately so. This makes them 
ideal for measuring the rotation rate of the lens. A feature 
occurs on the lightcurve when a magnification pattern fea- 
ture sweeps over the source. A radial feature that subtends 
the angle -i/jf relative to the binary axis will occur on the 
lightcurve when 



9{t)^<t>{t) + t^f. 



(15) 



By solving this equation we can use the timing of repeated 
features to easily obtain approximate measurements of some 
of the lens parameters. This means that many of the lens 
parameters can be 'read-off' the lightcurve, and it is pos- 
sible to build an approximate model of the lens quickly, 
without complex analysis. For such estimations, the most 
important magnification map features are the magnification 
arms (shown in Figure [l] which extend from the central 
caustic to the secondary caustics), and a planetary demag- 
nification, a feature that is only present for lenses with small 
mass ratios, which is a region of demagnification relative to 
the single lens that lies between the secondary caustics, with 



10 M. T. Penny, E. Kerins and S. Mao 




Figure 9. An example lightcurve by an RRL showing how 
hghtcurve features relate to the parameters of the lens. The lens 
has parameters tg = 30 d, T = 38 d, s = 0.3, q = 0.1, uq = 0.8, 
00 = 2.14 



its minimum lying along the binary axis. Both features are 
complementary, as in equal mass ratio binaries the plane- 
tary demagnification does not occur, but the magnification 
arms are strong and very close to radial, while in low mass 
ratio binaries, the magnification arms are weaker and less 
radial, but the demagnification region is strong. 

Figure [^ shows a lightcurve where features repeat 
strongly five times. The first step to estimating RRL pa- 
rameters is to fit the lightcurve with a Paczyiiski curve. This 
is relatively trivial, and most RRL lightcurves will approxi- 
mate a Paczyiiski curve with only small deviations. This fit 
allows an accurate estimation of the parameters to, ^e and 
|iio|, the last down to an ambiguity in sign, which corre- 
sponds to the source passing the lens centre on its left (pos- 
itive Mo) or right (negative uq) and imposing that the lens 
always rotates anti-clockwise. This Paczyhski model com- 
pletely describes the source trajectory, and hence defines the 
left hand side of equation [TS] The orbital period can now be 
estimated by timing two occurrences of the same magni- 
fication pattern feature. The period is not simply the time 
elapsed between two features, because the source moves dur- 
ing this time. Instead, by solving equation[T5]we can find the 
relation between the period T, and the time of two consec- 
utive occurrences of the same magnification pattern feature 
at times t\ and t2 

27t 

^- 2.+ [.(t.)-.(tOl (^-^-)' ^''^ 

where the fraction is the number of orbits the source com- 
pletes between the two source encounters. The degeneracy 
in the measurement of the sign of uq affects this equation, 
due to the presence of the 6{t) terms, but can be resolved 
if more than one pair of features is available for estimating 
T, as only one value of mq will give consistent estimates of 
T for different feature pairs. 

With an estimate of the period, if we know the angle 
subtended by a feature on the magnification map i/)f , we can 
also estimate the phase angle 0o, again taking into account 
the source motion 



27r 

ho = 0{tf)-ipi-—(t{-to) 



(17) 



makes this task simple. However, the demagnification may 
not be obvious, or if the mass ratio of the lens is high, may 
not be present. In these cases it is necessary to know ijjt 
for the magnification arms. Knowing that they extend from 
the central caustic (roughly at the centre of mass) to the 
secondary caustics, we need only know the po sition of the 
secondary caustics to estimate f/if . iBozzal l|2000l ) has derived 
analytical approximations to the position and shape of sec- 
ondary caustics in close lenses using a series expansion of 
the Jacobian, critical curves and caustics for s ^ 1. He 
finds that the secondary caustics are located at 

s(i + g) [ ±yg(2-s2) 

in the static centre of mass system. Figure [8] shows that 
this expression is reasonable even when s is quite large. If 
we assume the magnification arms are radial, we can use 
Equation [TS] to approximate the angle of the magnification 
arms, to second order in s, as 



(18) 



i/'i 



arctan 



±^(2 + g' 



(19) 



which is relatively insensitive to the lens separation s. It is 
useful to note the asymptotic behaviour: ip± ~ ±2g^" as 
g — >■ and ^})± — )■ ±7t as g — )■ 1. While the dependence of 
V'i on q implies an ambiguity in the estimation of 0o, the 
corollary is that we can estimate the mass ratio from the 
timing of features as well. Using the times of consecutive 
magnification arm crossings, t+ and t_, we have 



m - \ 



Kt-)~o{t+) 



27t , 

T ^ 



t+) 



(20) 



This value can be substituted into equation 1171 and equa- 
tion [19] can then be solved for g. 

The remaining parameter that we are interested in is 
the lens separation s. The angle of features is essentially in- 
dependent of s, so it is not possible to estimate s by timing 
features. However, by noting that the magnification pattern 
becomes essentially featureless beyond the secondary caus- 
tics (see Figure [l|, and that the position of the caustics does 
depend on s, it is possible to estimate s from the lightcurve. 
Unfortunately the secondary caustics are very small, and 
in most events, they will not pass directly over the source, 
so the estimate will not be very accurate. The best esti- 
mate of the position of the caustic will be derived from the 
largest peak due to a magnification arm in the wings of the 
lightcurve (e.g. the peak at i ~ —90 d in Figure|9]). This will 
occur when the radial source position approximately coin- 
cides with the radial caustic position, so that \zb\^ « l-^±P- 
Using Equation [6] to first order, we can write 

, . '1/2 
2 , I tc ~ to 



Wo + 



fE 



(21) 



The planetary demagnification region has tpt = 0, which 



where tc is the time of the peak due to the caustic. 

We have outlined how the parameters of an RRL can 
be estimated from pairs of feature timings in the case of the 
simplest RRL. However, in a given event there may be many 
repetitions, and better parameter estimates can be obtained 
by considering all the lightcurve features simultaneously. For 
a given magnification pattern and source trajectory it is pos- 
sible to compute a timing model by finding all possible so- 
lutions of equation 1151 6{t) = ^(t) -I- ?/>f for each feature. By 
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piece of information is needed to break the degeneracy. This 
can be done bymeasurine tte ~ au(l — x)/rE, the microlens- 
ing parallax (j Gould] [l993), or by measuring 6e ~ te/Di, 
the angular Ei nstein r adius , through detection of finite- 
source effects IIGould 1994 ; iNemiroff fc Wickramasinghd 
1 1994 IWitt fc Maolll994l '). or direct de tection of the lens 
once it has separated from the source (lAlcock et al.1 I2001I : 



_- — je pa — — — — — , - — , 

iBennett et allbood l2007l : iKozlowski et al-lbOO T''). Measure- 
ment of either tte or 9 e allow a seco nd mass-distance relation 
to be written, for tte l|Gouldlll993 ') 



M 



au^(l — x) 

C^xDsTvl ' 



or similarly for 6e 

INemiroff fc Wickramasinghelll99l ) 



l|Gouldl 



(23) 
[l99i; 



Figure 10. Plot of the various mass-distance relations for the 
event shown in Figure [9] labelled by the parameter measurement 
that would allow their definition. The arrows point into the region 
that is allowed should only an upper limit on T, tte or 6e be 
available. If the period T is measured along with only one of tte 
or 9e, the mass and distance to the lens can not be determined 
uniquely, but even a relatively weak upper limit on the other 
parameter may be sufficient to rule out one possible solution; 
note however that a lack of finite-source effects places a lower 
limit on 6e- 



extracting the occurrence time of all the lightcurve features 
it is possible to fit timing models to this timing data. It is 
also possible to add additional features to this timing model, 
such as the effects of inclination and eccentricity by modify- 
ing the function <^(f), or microlensing parallax by modifying 
9{t). This modelling may be significantly faster than a full 
lightcurve fitting analysis, especially when additional effects 
are included, as there is no need to calculate finite-source 
magnifications. While it will not fully remove the need for 
lightcurve fitting, it will significantly narrow down the range 
of parameters over which lightcurve fitting has to search. 

We have shown that it is possible to estimate the pa- 
rameters of an RRL lightcurve, but what we would really 
like is to be able to measure the physical parameters of the 
lens, most importantly the lens mass and the binary separa- 
tion in physical units. Compared to a static binary lens, we 
have one additional piece of information with which to infer 
M and a: the orbital period. iDominikj (|1998i ) has shown that 
by combining the orbital period T and the lens separation 
s, it is possible to write down a mass-distance relation 



M = 



C2(l 



(24) 



M = 






C6s6x3(l-X)373r 



(22) 



which relates the mass to the lens distance through 
known quantities, assuming the source distance is known 
from its colour and magnitude; the constant C — 
2.85M^ au kpc~^ " when the period is measured in years 
and the source distance in kpc. As demonstrated in Ap- 
pendix IB3I it is likely that if the orbit is inclined it will be 
possible to measure the inclination and account for projec- 
tion, so that the value of s measured is a good approximation 
of a/vE ■ This means that as equation [22] has a minimum at 
X — 0.5, we can place a firm lower limit on the mass of the 
lens, and an upper limit on the semi-major axis. 

To improve on the mass-distance relation, an additional 



if 6e is measured in mas. One of these relations can then be 
combined with equation[21]to yield two possible solutions to 
the mass and distance. This can be seen in Figure [Tol which 
plots the mass-distance relations for the event shown in Fig- 
ure J^J The tte- and 6'E-lines cross the T-line in two places: 
once at the true parameter values x = 0.95, M — IMq, 
and once at other values of M and x which are different 
for each relation. With a measurement of only one of tte or 
9e, it is not possible to uniquely determine the mass and 
the distance. This is likely to be the case, as finite-source 
effects are most likely in lenses close to the source, while 
parallax is most likely in lenses close to the observer. How- 
ever, even a crude limit on the unmeasured parameter may 
be enough to rule out one possible solution; e.g. an upper 
limit on tte from the lack of parallax effects may allow the 
solution with smaller x to be ruled out, or a lower limit on 
9e from the lack of finite-source effects may allow the solu- 
tion with larger x to be ruled out. Direct detection of the 
lens may require a very long time baseline as RRL features 
are most detectable in events with low lens-source proper 
motions. However, RRLs are more likely to be more massive 
than the average lens, and hence brighter, and the diffraction 
limit of 30-40 m class telescopes, such as the Thirty Metre 
Telescope (TMT)|3the Giant Magellan Telescope (GMT)|3 
or the European Extremely Large Telescope (E-ELT) ]j may 
be sufficient to resolve the lens and source in a reasonable 
time. 



5 DISCUSSION AND CONCLUSION 

Although the phenomena of microlensing by lenses with 
rapid orbit al motion has been discussed in the literature 
previously l|Dominik]] 19981 : iDubath et al.ll2007l '). the realities 
required to estimate its occurrence rate have not properly 
been accounted for. In this paper, we have established a 
firm theoretical basis for rapidly rotating lenses, and used 
it to estimate the range of parameters over which they are 



^ Note that parallax or finite-source effects were not included in 

the model used to plot the lightcurve 

2 http://www.tmt.org 

^ http://www.gmto.org 

^ http://www.eso.org/public/teles-instr/e-elt.html 
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detectable and the rate that they are expected to be ob- 
served. We find that RRLs with masses and orbital radii 
typical of binary stars are detectable, and that there is a 
reasonable chance that they will be detected, either in cur- 
rent microlensing data sets or in ongoing or near-future mi- 
crolensing surveys. 

In calculating these rates we have actually used the 
relatively stringent criteria of requiring that two or more 
lightcurve features from the same orbital phase are detected 
in the lightcurve. If we relax this repetition requirement 
somewhat, to include lenses that display significant signs 
of orbital motion (say several degrees rotation per is), the 
event rate will increase significantly, as lenses can then have 
larger orbits and hence stronger lightcurve features. We have 
shown in a previous paper that orbital motion is detectable 
in a large fraction (~ 15 percent) of binary lenses with de- 
tectable binary lensing features and orbital periods co mpa- 
rable to the microlensing timescale IjPennv et al.ll201H ). 

We have detailed how the features of an RRL lightcurve 
can be used to measure its period and potentially measure 
its mass. Even if features do not repeat, if several features 
are detectable in the lightcurves of binary lens events, then 
the techniques we have outlined for timing features and ex- 
tracting parameter estimates may be of some use in their 
analysis. Without repeating features, the orbital period may 
not be constrained as accurately, but it should be possible to 
place constraints on the lens mass and orbit in many cases. 

So far, we have neglected to discuss the prospects for 
positively identifying RRL events from other events which 
may mimic their features. Periodic features may also be in- 
duced by orbital motion in the observer and source planes, 
or intrinsic variability in the source or a blend star. In the 
observer plane, the period of orbital parallax effects is well 
defined, and unless the lens has an orbital period similar 
to 1 yr it is unlikely the effects will be confused. Even if 
the orbital period is close to one year, the shape of features 
in the lightcurve are likely to be different. Orbital effects 
in the source plane may be more difficult to exclude, as 
the period is not fixed. If t here is only a s i ngle luminous 
source (the xallarap case. |Paczvnskil 119971 : iHan fc Gouldl 
ll997l : lRahvar fc Dominikll2009l ). a timing analysis similar to 
the one we proposed for the lens can be performed for the 
source. This analysis is somewhat easier and more precise for 
xallarap as there are no complicated features in the magni- 
fication pattern. If this analysis is insufficient to separate 
the two cases then the shape of lightcurve features should 
differentiate the two interpretations. In the case where both 
sources are luminous, the lightcurve can take a more com- 
plicated s hape, which may more close l y resemble that of a n 
RRL (e.g. ICherepashchuk et ai]|l995l : JHan fc Gouldl 1 199'il ). 
Even in this case, timing analysis for maxima and minima 
of the lightcurve should be easier than for RRLs, and full 
lightcurve modelling starting from timing analysis solutions 
will likely be able to differentiate the two scenarios. Finally, 
source variability that is not detectable at baseline, but be- 
comes detectable with the increased photometric accuracy 
thanks to microlensing magnification may also produce sim- 
ilar lightcurve features. 

It is worth noting that we should naively expect the 
rate of RRL/significant lens orbital motion events to be 
similar to the rate of binary source orbital motion events, 
as the factors that govern their occurrence, such as the ra- 



tio of orbital separation to the Einstein ring, and the ratio 
of orbital to microlensing timescales will have similar dis- 
tributions in the lens and source populations. Naively, we 
would expect the rate of parallax events to be roughly ten 
times greater than the rate of RRL events with orbital pe- 
riods ~ 1 yr, as the binary fraction is ~ 0.1 while the ob- 
server is always orbiting. It is worth comparing this with 
the nu mber of repo rted single lens paralla x events, ~ 20- 
50 (e.g. |Poindexter~et al. 20051 : ISmith et aL|j2005 . and refer- 
ences therein), while ~ 10 events have been successfully fit 
with xallarap models IjSmith et al.l 120031 : IPoindexter et al.l 
[2OO5I). In contrast, only one binary lens event has been suc- 
cessfully fit with a binary lens model that shows sig nificant 
rotation, MACHO-97-BLG-41 (JAlbrow et all I2OOOI ), which 
rotates at a rate ~ 4° per tE, a low rate compared to RRLs, 
that is detected thanks to the source crossing the central 
and one secondary caustic, as opposed to the smaller, more 
continuous features of RRLs. It is possible therefore that 
many events with significant rotational orbital motion sig- 
natures have not been modelled, or have been interpreted 
as xallarap events. It is thus important that any event that 
is modelled with xallarap also be tested with an orbiting 
binary lens model. 
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APPENDIX A: WHAT IS HAPPENING TO THE 
IMAGES? 

The image configuration of a point-mass lens consists of two 
images: a major image, of positive parity and magnification 
/^+ > 1, outside the Einstein ring and a minor image of 
negative parit y and magnification ^, < 0, inside the Ein- 
stein ring (e.g. lRefsdailll964l : lLiebeslll964l '). The addition of 
a second mass to the lens causes an additional image of neg- 
ative parity t o be produc ed if the source does not lie within 
a caustic (S chneider fc W eiss 1986). If the lens is far from 
resonance, i.e. s ^ 1 or s 3> 1, two of the three images can 
still be associated with the major and minor images of the 
single lens, while the new third image is labelled a tertiary 
image. 

It is interesting to study what is happening to each 
of the three images during t he co urse of an RRL event. 
iDubath. Gasparini fc Durrerl (|2007|) study the effects of an 
orbiting close binary lens on the major image, by casting 
the lensing potential as a time varying quadrupole. They 
show that the major image can exhibit significant time de- 
pendent deviations from the single lens form when it is 
highly magnified, and go on to calculate the expected rate 
of events showing such deviations. Unfortunately, they ne- 
glect to treat both the tertiary image and the minor image, 
the latter of which will be magnified to a simil ar degree, as 
Im+I - Im-I = 1 (e.g lRefsdall[l964l : lLiebeslll964l l. 

In Figure lAll we plot the lightcurves of all three im- 
ages for an unrealistic lens with repeating features. In the 
top panel of the figure the combined lightcurve (that which 
would be observed) is shown in red, and the absolute mag- 
nification of the major j^+|, nunor |^_| and tertiary im- 
ages 1/^3 1 is plotted in green, blue and magenta respectively. 
The second panel from the top shows the residual of the 
lightcurve relative to the lightcurve of a single lens of the 
same total mass, and the third and fourth panels show the 
absolute magnification residuals of the major and minor im- 
ages with respect to their single lens counterparts; the bot- 
tom panel shows the absolute magnification of the tertiary 
image (note the different scale) , which does not have a coun- 
terpart for the single lens. Each image shows a strikingly 
different pattern of features: the major image is only signif- 
icantly perturbed from its single lens form when the source 
is within ~ te of the centre of mass, while the minor image 
shows significant perturbations out to the position of the 
secondary caustics. It is only when the source is close to the 
secondary caustics that the tertiary image is magnified sig- 
nifican tly. However, the mo st important aspects in relation 
to the lDubath et al.l (|2007h event rate calculation are that 
firstly, the minor image can experience larger perturbations 
than the primary image, and secondly that the perturba- 
tions of the major and minor images are out of phase by 
7t, in such a way that the peak to trough amplitude of the 
deviations in the total amplification lightcurve can be sig- 
nificantly reduced from those of the lightcurves of individual 
images. 



APPENDIX B: ADDITIONAL FACTORS 
AFFECTING RRL DETECTABILITY 

In the main text we have mentioned a number of addi- 
tional effects that can affect the form of an RRL lightcurve 
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and its detectability. Here we briefly outline the three most 
important effects and the impact they will have on RRL 
lightcurves and detectability. 



Bl Blending 

For a given photometric precision Gm magnitudes, the effec- 
tive threshold at the event baseline is 



Ammin 



2.5 log (10" 



1 + /,) -2.5 log /s 



2.5 log 1 + 0.922^ 



(Bl) 



where the approximation applies for small cjm and /^ is the 
fraction of the total light at baseline contributed by the un- 
lensed source. Figure iBll shows this for various values of the 
photometric threshold. It is clear that only with the most 
accurate photometry will it be possible to detect RRL fea- 
tures when the blend contributes most of the flux, and for 
less accurate photometry, am ~ 0.02 even a small amount 
of blending will significantly affect the detectability of fea- 
tures. The effect of blending decreases as the magnifica- 
tion increases, but we wish to see features over the entire 
lightcurve, and only a small region of the lightcurve will be 
magnified enough to significantly reduce the effect of blend- 
ing. 



B2 Finite-source effects 

Figure IB2I shows the lightcurve of an RRL lensing a giant 
source of radius 100-Rq, in comparison to the same RRL 
lensing a point source. The effect of the finite source on the 
lightcurve is clear, as it causes a wider, lower peak magnifica- 
tion. Whilst the lens centre of mass transits the source, there 
is effectively no deviation from the finite-source point-lens 
lightcurve, except for spikes in the residual at f ~ ±20 d, 
which are char acteristic of a large source cr ossing a small 
central caustic (|Dong et al.ll2009l : lHanll2009l 'l. In the wings 
of the lightcurve there is very little difference between the 
finite- and point-source lightcurves, and most of the fea- 
tures in the residuals have the same amplitude. Only when 
the source is very close to the secondary caustic is there 
any deviation from the point-source lightcurve in the wings. 
The left inset of Figure [B2l shows that the peak in the finite- 
source lightcurve at f « —210 d is slightly broader and about 
half the amplitude of the point-source lightcurve. Interest- 
ingly, this peak, although broadened by the finite-source, is 
still much narrower than the source crossing time, which de- 
termines the width of the central peak. Its width is instead 
determined by the time taken for the secondary caustic to 
cross a source diameter. 

The example we have shown is very extreme, with a very 
large source, very close to the lens, and even then the finite- 
source effects only render binary features undetectable over 
a relatively small fraction of the lightcurve. A typical giant 
source star will be up to a factor of ten smaller, so the part of 
the lightcurve severely affected by finite-source effects will be 
correspondingly smaller. As the source has to be transited by 
the lens centre for finite-source effects to become apparent 
at the lightcurve peak, the probability of this occurring is 
also reduced by the same factor. This means that finite- 
source effects will not affect the detectability of repeating 
features very much. If finite-source effects are detected in an 



event, the measurement of the source radius, combined with 
a measurement of the lens period can be used to gether to 
measure the lens mass to a two fold degeneracy (jDominilg 
[19981 ). 



B3 Inclination and eccentricity 

Inclination and eccentricity of the lens orbit will act to make 
the magnification pattern motion much more complicated, 
as changes in the projected lens separation cause the caus- 
tics to move and change shape (see e.g. Figures 21 and 22 of 
Penny et al. 2011 ). The effects are too complicated to inves- 
tigate in detail, but it is worth considering them in brief. For 
a lens with a given semimajor axis, inclining the orbit should 
reduce the detectability of features over part of the orbit, as 
s decreases. Figure [B3l shows the effect of inclination on the 
lightcurve of an RRL. It shows that inclination tends to de- 
crease the amplitude of features, but does not completely 
wipe them out, even when the inclination i = 90°. In this 
extreme case, rather than rotating, the secondary caustics 
move along diagonal lines as the projected separation of the 
lenses changes, but their angle does not, except for fiips by n 
every half period. Inclination significantly changes the mor- 
phology of the lightcurve, and can also change the timing of 
peaks (see e.g. those at t ~ —300 d), which implies that it 
may be possible to measure the inclination of the lens orbit 
from the lightcurve. 

In contrast to inclination, eccentricity may increase the 
detectability of features. For a given semimajor axis, eccen- 
tricity can both increase and decrease the projected sepa- 
ration. However, Kepler's second law implies that the lens 
will spend longer at larger projected separations (assuming 
no inclination). As with inclination, eccentricity will also 
change the lightcurve morphology and timing of features, so 
it may also be possible to measure the eccentricity of the lens 
from the lightcurve. Simultaneously including the effects of 
inclination and eccentricity in the modelling of an RRL event 
will likely be difficult, as together they require an additional 
four parameters over the standard RRL parametrization. 
However, as the angle of magnification pattern features does 
not depend strongly on the projected separation, it will be 
possible to include inclination and eccentricity in the tim- 
ing analysis proposed in Section |4l which may significantly 
ease the analysis by narrowing down the search space to the 
range of parameters compatible with timing measurements. 



B4 Parallax 

Parallax effects due to the motion of the Earth about the 
Sun will cause the source to appear to take a curved path 
through the plane of the sky, and will affect the lightcurve 
of an RRL event. If the magnitude of the parallax effect 
is small then it will cause only small perturbations to the 
shape of the lightcurve and the timing of features. Larger 
effects may cause significant changes to the RRL lightcurve, 
significantly changing the timing of features, and possibly 
making them appear less periodic, or adding a stronger an- 
nual periodicity to the lightcurve. However, while parallax 
may significantly complicate the interpretation of an RRL 
event, it does not affect the magnification map, and the de- 
tectability of RRL features should remain the same. More- 
over, the detection of parallax in an RRL event will allow the 
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lens mass to be measured, at least to a two-fold degeneracy 
(see Section |4]). Due to the photometric accuracy required 
to detect RRLs, and the long timescales of the events, the 
probability of detecting parallax along with R RL features is 
significant (jBuchalter fc Kamionkowskilll997^ . 
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Figure Al. Liglitcurves and residuals for each image of a mi- 
crolensing event with repeating features. The top panel shows 
the absolute magnification of the combined images (|/i|), and the 
individual major (|/^+|), minor {|/i-|) and tertiary ([^sD images 
in different shades of grey. The central panels show the absolute 
magnification residual with respect to the single lens form for all 
images combined, the major image and the minor image, going 
from top to bottom, respectively; the bottom panel shows the 
absolute magnification of the tertiary image, which has no sin- 
gle lens counterpart. The event has the unrealistic parameters 
no = 0.4, s = 0.3, q = 1.0 and ts/T = 10. 
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Figure Bl. The effect of blending on the photometric detection 
threshold. The effective threshold Amniin is plotted against the 
ratio of source to total blend flux /a for three values of photomet- 
ric precision am- The solid lines show the exact value, whereas 
the dashed line shows the approximation for small am- 




-100 

Time (days) 



300 



Figure B2. The lightcurve of an RRL lonsing a finite source 
of radius IOO-Rq (black) compared to the lightcurve of the same 
RRL lensing a point source (grey). The inset figures show in more 
detail the residuals when the source is close to the secondary caus- 
tic (on the left) and the central caustic (on the right). The lens 
has a mass M = O.SMq, semimajor axis a = 0.4 au, mass ra- 
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Dg = 8 kpc, source velocity Vt = 50 km s~^, impact parameter 
uo = 0.1 and phase angle 0o = 7r/4. The ratio of source to Ein- 
stein angular radii pg = 9^/0^ = 0.28 is very large. The effects 
of finite sources are only significant when the source is near the 
central or secondary caustics. 
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Figure B3. The lightcurves of RRLs with different orbital in- 
clinations relative to the line of sight. For each lightcurve, the 
lens has mass M = 0.58Mq, semimajor axis a = 0.54 au, mass 
ratio q = 0.52, fractional lens distance x = 0.86, source distance 
Db = 9.5 kpc, source velocity vt = 61 km s~^, impact parameter 
MO = 0.77, phase angle (f>o = 4.3 measured in the plane of the 
orbit. The orbit was circular, and inclined about the x'-axis as 
defined in Figure [S] 



